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ABSTRACT 


The  results  of  a  series  of  field  trials  on  the  diffusion  and 
ground  deposition  of  100  micron  glass  microspheres  from  a  continuous  point 
source  at  a  height  of  92  metres  are  discussed.  The  observed  crosswind 
integrated  deposit  density  as  a  function  of  distance  from  the  source  was 
used  to  test  two  prediction  models.  One  of  these  models  employs  appro¬ 
priately  averaged  standard  deviations  of  vertical  turbulence  as  the  main 
parameter  of  atmospheric  diffusion.  The  other  is  the  steady  state  K- 
Theory  diffusion  model  with  a  coefficient  of  eddy  diffusivity  which  varies 
with  height.  In  general,  there  was  reasonably  good  agreement  between  the 
observed  and  predicted  crosswind  integrated  deposit  density  as  a  function 
of  distance,  for  the  sloping  plume  model.  However,  the  K-Theory  model 
predicts  a  peak  deposit  much  lower  than  observed  and  a  more  gradual  decrease 
in  the  deposit  density  than  observed. 
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DIFFUSION  AND  GROUND  DEPOSITION  OF  100  MICRON  PARTICLES  FROM 


by 

0.  Johnson 

INTRODUCTION 

During  the  period  between  1958  and  1967  several  series  of  field 
trials  were  conducted  at  the  Defence  Research  Establishment  Suffield  in 
a  study  of  the  diffusion  and  deposition  on  the  ground  of  particles  released 
from  low  level  sources.  Hage  { 1961 )  reported  on  a  series  of  field  trials 
in  which  glass  microspheres  of  100  microns  diameter  were  released  from 
a  continuous  point  source  at  a  height  of  15  metres,  and  Walker  (1965) 
reported  on  similar  trials  using  50  micron  particles  released  from  8  metres. 
Results  of  seven  field  trials  in  which  30  micron  glass  microspheres  were 
released  from  a  height  of  2.75  metres  have  been  reported  by  Johnson,  McCallum 
and  Larson  (197**). 

For  these  series  of  field  trials,  mathematical  models  for  predicting 
the  downwind  distribution  of  deposit  density  on  the  ground  have  been  tested. 

In  addition  the  mass  recovery,  estimated  from  the  deposits  on  the  surface 
samplers,  was  compared  with  the  total  mass  emitted. 

In  this  paper  the  results  of  lU  field  trials  on  the  diffusion 
and  deposition  on  the  ground  of  glass  microspheres  of  nominal  100  micron 
mass  median  diameter  from  a  continuous  point  source  at  a  height  of  92  metres 
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are  discussed.  (Trial  2  was  unsuccessful.)  The  main  objective  of  this  study 
was  to  provide  a  mathematical  model  which  adequately  predicts  the  downwind 
distribution  of  the  deposit  density  of  particles  of  terminal  velocities  equi¬ 
valent  to  those  of  about  100  micron  glass  microspheres.  Two  models  were 
tested:  a  Gaussian  Sloping  Plume  model.  Walker  (1965),  and  a  Steady  State 
K-Theory  model,  Monaghan  and  McPherson  (private  communication). 

EXPERIMENTAL  DETAILS 

The  sampling  technique  and  meteorological  instrumentation  were 
similar  to  those  described  by  Hage,  Diehl  and  Dudley  (i960).  Glass  micro¬ 
spheres  coated  with  fluorescent  dye  were  emitted  continuously  for  periods 
of  28  to  60  minutes  from  the  top  of  a  92  metre  tower.  The  particles  were 
collected  on  horizontal  adhesive  sampling  surfaces  placed  on  the  ground 
along  crosswind  arcs  at  downwind  distances  up  to  It, 828  metres.  The  number 
of  particles  collected  on  each  sampling  surface  were  counted  visually  under 
ultra-violet  light,  by  means  of  a  magnifier. 

The  distribution  of  particle  sizes  by  mass  was  determined  for  two 
samples  taken  from  the  lot  used  in  the  field  trials.  The  data  are  given 
in  Table  1. 

On  four  of  the  field  trials,  microspheres  from  the  sampling  surfaces 
on  a  small  number  of  arcs  were  measured  by  microscope,  so  that  the  variation 
in  mass  median  diameter  with  downwind  distance  could  be  compared  with  that 
predicted  from  the  model,  using  the  distribution  in  the  bulk  sample.  The 
data  are  given  in  Table  2. 

Meteorological  measurements  in  support  of  the  field  trials  were 
made  near  the  source  for  a  period  of  60  minutes,  which  included  the  period 
of  emission  (Table  3).  Wind  speedB  were  recorded  at  8  levels  on  the  92  metre 
tower  (Table  h)  and  the  temperature  difference  between  heights  of  13  and  1  metre 
were  measured  (Table  5)-  Crosswind  and  vertical  angular  deflections  of  light 
vanes,  mounted  at  heights  of  16,  1*8  and  92  metres  on  the  tower,  were  recorded 
continuously  on  magnetic  tape  and  high  speed  recorder  paper  charts.  The 
magnetic  tape  record  was  digitized  at  one  second  intervals  and  the  standard 
deviations  of  the  vane  angles  were  computed  with  data  averaging  times  of 
1,  5,  10,  20,  60  and  100  seconds. 
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PREDICTION  MODELS  FOR  DOWNWIND  DISTRIBUTION  OF  DEPOSIT  DENSITY 
(a)  Sloping  Plume  Model 

Since  the  microspheres  used  for  these  field  trials  vary  somewhat 
in  size,  there  is  a  certain  amount  of  downwind  spread  in  the  deposit  pattern 
which  has  to  be  considered.  Walker  (1965)  developed  a  prediction  model  in 
which  particles  in  each  size  class  were  spread  about  the  non-turbulent 
gravity  fall  path  by  Gaussian  vertical  turbulence,  with  angular  standard 
deviation  a.  The  horizontal  wind  speed  IT  was  considered  constant,  and  if 
the  terminal  velocity  of  the  particle  size  class  is  'v* ,  then  the  distri¬ 
bution  of  mass  Q  about  the  gravity  fall  path  is: 

=  — 9 —  exp  (_e2/2o2) .  (l) 

de 

where  9  =»  <p  -  tan-1  — 

u 

and  0  is  the  angular  standard  deviation  of  Gaussian  turbulent  spread. 

These  relationships  are  shown  in  Figure  29.  The  distance  downwind  'x' 
travelled  by  a  particle  before  landing  is 

x  =  h  cot  $ 

Hence  0  *  tan-1  —  -  tan  -1  ~  ~  ~ 

x  u  x  TT 

The  deposit  density  is  then 

»  -  If  I  ■  ©  ®  *  p  7^=r  MW> -  (*> 

*  v  d7T  o 
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The  average  of  the  standard  deviation  o  at  l6,  U8  and  92  metres 
was  taken  from  Table  6,  as  equivalent  to  the  vertical  vane  angle  standard 
deviation  when  the  data  were  averaged  over  periods  equal  to  the  travel  time 
divided  by  the  factor  g,  which  represents  the  ratio  between  Lagrangian  and 
Eulerian  time  scales,  Pasquill  (1962).  The  travel  time  vas  taken  as  the 
distance  to  each  sampling  arc  divided  by  the  mean  of  the  wind  speed  at 
l6,  U8  and  92  metres.  8=1  and  8  =  were  used  in  the  prediction  model. 

(b)  Steady  State  K-Theory  Model 

In  the  K-theory  model,  Monaghan  and  McPherson  (private  com¬ 
munication),  the  atmospheric  dispersion  of  particulate  material  is  described 
by  the  two-dimensional  steady  state  diffusion  equation 


u(z)  3c/3x  =  3{K(z)  3c/3z)  /3z  +  q3c/3z 


(3) 


where 


c  =  steady  state  concentration  at  x,z  for  a  continuous 
line  source,  or  total  dosage  when  the  line  source 
is  instantaneous . 

K(z)  =  the  coefficient  of  eddy  diffusivity  at  height  z. 

u(z)  =  the  mean  horizontal  wind  speed  at  height  z. 

q  =  terminal  velocity  of  the  particles,  assumed  monodisperse. 


It  is  assumed  that  there  is  an  impervious  "lid"  for  diffusion  at  height  H, 
where  the  vertical  flux  of  material  is  zero  or, 


Lim  F  =  0 
z-»H 

The  flux  of  material  at  the  lower  icindary  of  the  turbulent  atmosphere  is 
given  by 


Lim  {-F}  =  E(P&  +  q/S)S  Lim  c 
z-O  z-K) 
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Pr  is  the  mean  transport  velocity  of  non-settling  material  between  the 
rough  surface  and  the  turbulent  boundary  layer  and  is  related  to  the 
sublayer  Stanton  member  B  by 

P  =  u  B/S 

a  * 

where  S  is  the  total  specific  surface  of  the  roughness  elements,  and  E 
is  the  efficiency  of  retention  of  the  surface. 

The  boundary  conditions  are  therefore 

Lim  (K(z)  3c/3z  +  qc}  =  0 
Z-+H 

and  Lim  (K(z)  3c/3 r  +  qc}  =  E(Pa  +  q/S)S  Lim  c, 
z-*0  z-*0 

with  the  initial  condition 

Lim  c  =  Q6 ( z-h)/u(h) 
x->-0 

where  Q  is  the  source  strength,  h  is  the  release  height  and  6(z)  is 
the  Dirac  delta  function. 

Simple  functional  forms  for  mean  wind  speed  and  the  eddy 
diffusion  coefficient  were  adopted  as  follows: 

u(z)  *  u(2)  {(x+A£.)/2}P 

and  K(z)  =  A(z+A!,}  u(2)  ; 

p  and  A  are  constants  related  to  atmospheric  stability  and  Ai.  is  a  length 
of  the  same  order  as  the  roughness  length  zq  of  the  surface.  u(2)  is  the 
mean  wind  speed  at  2  meters. 

In  the  above  equations  the  following  assumptions  were  made: 

E  =  1 
S  =  1 

H  *  500  meters 
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P  =  0.01  u( 2) 
a 

A  =  Q.Oh 

Equation  (3)  is  solved  by  a  finite  difference  approximation  method 
in  which  the  discretization  in  the  vertical  is  in  term-:  of  equal  intervals  of 
diffusive  resistance  rather  than  length. 

DEPOSIT  DATA  AND  ANALYSIS 

The  particle  counts  on  each  sampling  surface  were  converted  to 
deposit  densities  in  mg/m2  by  using  a  mean  particle  mass  computed  from 
Table  1.  The  data  are  given  in  Appendix  A.  The  downwind  distribution  of 
deposit  density  was  determined  by  crosswind  integration  along  each  arc. 

This  was  done  by  multiplying  the  deposit  density  by  the  spacing  of  the 
samplers  and  then  summing.  Table  ^  shows  the  crosswind  integrated  deposit 
density  (CWID)  normalized  by  dividing  by  the  total  mass  emitted. 

The  CWID  at  each  sampling  distance  was  predicted  from  equation  (2) 
of  the  sloping  plume  model.  The  predicted  and  observed  normalized  CWID  as 
a  function  of  distance  for  6=1  and  6  =  1*  are  given  in  Table  7  as  Model  A 
and  plotted  in  Figures  1  to  1*+.  Predictions  of  the  CWID  distribution  were 
also  made  using  the  steady  state  K-Theory  model. 

Comparison  of  the  predicted  and  observed  distances  to  peak  deposit, 
peak  deposit  and  recovery  to  the  farthest  sampling  distance  are  given  in 
Table  8. 

DISCUSSION 


Table  8  shows  the  predicted  and  observed  deposit  density 
characteristics  for  all  the  trials.  In  general,  there  is  fairly  good 
agreement  between  predicted  and  observed  deposit  density  characteristics. 

The  average  of  the  predicted  distances  to  the  peak  deposit  was  9%  less  than 
observed  when  6  =  was  used.  The  average  of  the  predicted  peak  CWID  was 
19 %  higher  than  observed.  Predicted  recovery  to  the  maximum  distance  sampled 
was  15?  higher  than  observed.  The  reduced  observed  CWID  and  recovery  may  be 
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attributed  to  the  very  low  recoveries  for  trials  5,  7  and  8.  The  wind  speeds 
for  these  trials  were  very  high  and  may  have  raised  dust  which  subsequently 
settled  on  the  sampling  surfaces. 

Figures  1  to  lL  show  that  the  choice  of  6  =  1  or  3  =  U  makes  a 

considerable  difference  in  the  predicted  pattern  of  the  deposit  density  as 
a  function  of  distance  from  the  source,  particularly  with  respect  to  the 
distance  to  the  peak  deposit  and  the  magnitude  of  the  peak.  In  general 
the  observed  data  are  better  predicted  when  3  =  ^  is  used  in  the  sloping 
plume  model. 

The  results  of  the  sizing  of  the  particles  on  the  sampling 
surfaces  given  in  Table  2  show  that  there  is  only  a  small  decrease  in  the 
mass  median  diameter  with  distance  from  the  source.  In  view  of  these  results, 
the  observed  deposit  density  was  not  corrected  for  variation  in  particle  size 
with  distance  from  the  source. 

(b )  Steady  State  K-Theory  Model 

The  downwind  distribution  of  the  CWID  predicted  by  the  K-Theory 
Model  was  computed  from  equation  (3)  and  most  of  the  characteristics  of  the 
distribution  are  given  in  Table  8  as  Model  B.  The  predicted  and  observed 
characteristics  are  compared  as  well  as  the  predictions  made  by  the  two 
models . 

When  averaged  for  all  the  trials  the  ratio  of  the  predicted  to 
observed  distance  to  the  peak  deposit  was  0.9  for  the  sloping  plume  model 
(Model  A)  and  1.0  for  the  K-Theory  (Model  B).  The  ratio  of  predicted  to 
observed  peak  deposit  density  was  1.2  for  Model  A  and  0.^5  for  Model  B. 

The  ratio  of  predicted  to  observed  CWID  at  the  maximum  distance  sampled 
was  0.75  for  Model  A  and  2.0  for  Model  B.  The  ratio  of  predicted  to 
observed  recovery  to  the  maximum  distance  sampled  was  1.15  for  Model  A  and 
0.90  for  Model  B. 

The  distance  to  peak  deposit  and  recovery  are,  on  average,  well 
predicted  by  both  models.  However,  the  K-Theory  Model  predicts  a  peak 
deposit  much  lower  than  observed  and  a  more  gradual  decrease  in  the  deposit 
density  with  distance  than  observed. 
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CROSSWIND  SPREAD  OF  THE  DEPOSIT  DENSITY 

The  angular  standard  deviations  for  the  particle  deposit  density 
distribution  across  each  sampling  arc,  o^,  were  compared  with  the  crosswind 
vane  angle  standard  deviation,  o^,  using  data  averaging  times  equal  to 
the  travel  time  to  each  arc  (0  =  l)  and  one  quarter  of  the  travel  time 
(6  =  4).  The  data  are  given  in  Table  9  and  plotted  in  Figures  15  -  28. 

With  0  =  4  the  average  of  the  ratio  for  all  the  trials  was 


1.06.  For  9  of  the  trials  o  was  greater  than  o  .  _y_  varied  between  0.66 

y  p  o 

p 

and  1.35  except  Trial  7  for  which  the  ratio  was  1.49.  There  appears  to  be  no 

obvious  explanation  for  the  high  ratio  for  this  trial.  The  decrease  in  the 

observed  standard  deviation  with  distance  is  consistent  with  predictions. 

Similar  results  were  reported  by  Walker  (1965)  for  50  micron 

microspheres  released  from  a  height  of  8  metres .  For  30  micron  microspheres 

emitted  from  a  height  of  2.75  metres,  Johnson  et  al.  (1974),  the  ratio  0  /o 

y  P 

varied  between  0.73  and  0.83. 

CONCLUSIONS 

The  observed  crosswind  integrated  deposit  density  as  a  function 
of  distance  was  compared  with  that  predicted  by  the  Gaussian  sloping  plume 
model,  using  standard  deviations  of  the  plume  equivalent  to  the  vertical 
turbulence  standard  deviations  when  the  data  were  averaged  over  periods 
equivalent  to  the  travel  time  to  each  sampling  arc  (0  =  l)  and  one  quarter 
of  the  travel  time  (0  =  4).  For  most  of  the  trials  there  was  good  agreement 
between  the  observed  and  predicted  deposit  density  distribution  when  0=4 
was  used. 

The  distance  to  the  peak  deposit  and  recovery  were  well  predicted 
by  the  Steady  State  K-Theory  model.  However,  the  predicted  peak  deposit 
was  much  lower  than  observed  and  there  was  a  more  gradual  decrease  in 
the  deposit  density  with  distance  than  observed. 


UNCLASSIFIED 


V 


* 


UNCLASSIFIED 


/9 


A  comparison  of  predictions  by  the  two  models  indicates  that 
the  data  are  better  predicted  by  the  sloping  model  than  by  the  K-Theory 
model. 

The  standard  deviation  of  the  crosswind  distribution  of  the 
deposit  density  was  compared  with  the  standard  deviation  of  crosswind 
turbulence  using  averaging  periods  for  0=1  and  0  =  1+.  There  was 
reasonably  good  agreement  between  observed  and  predicted  standard  deviations 
when  0  =  1+  was  used. 
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